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Abstract 

An  existing  method  of  calculating  the  elgenmodes  of  confocal  unstable 
strip  resonators  is  reviewed.  This  method  has  the  advantage  of  accurately 
accounting  for  diffraction  effects  but  is  limited  to  inclusion  of  an 
arbitrary  symmetric  gain  distribution.  An  extension  of  this  method  to 
Include  a  realistic  saturated  gain  distribution  is  proposed.  This  gain 
distribution  is  calculated  numerically  by  using  an  Iterative  technique  to 
arrive  at  a  self-consistent  solution  for  both  the  gain  and  field  distribu¬ 
tions.  Several  appendicles  are  presented  to  aid  In  understanding  the 
basic  theory. 
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STUDY  OF  APPLICABILITY  OF  ITERATIVE  SATURATED 
GAIN  MODEL  IN  CONFOCAL  UNSTABLE  RESONATOR  ANALYSIS 


I.  Introduction 


Background 

Unstable  laser  resonators  have  long  been  recognized  as  desirable 
cavities  for  use  in  high  power  lasers  because  of  their  large  mode  volume, 
direct  output  coupling,  and  mode  control  properties  (Ref  8).  As  the  use 
of  unstable  resonators  has  expanded,  the  need  to  develop  quick  and  accu¬ 
rate  methods  of  analytically  studying  the  behavior  of  the  modes  within 
the  resonator  has  become  more  important.  Since  the  positive  branch  con- 
focal  resonator  is  an  important  resonator  for  many  applications  and 
representative  of  this  class  of  resonator,  its  analysis  is  considered  in 
this  paper. 

The  primary  difficulties  in  analyzing  unstable  resonators  arise  from 
the  rapidly  oscillating  field  caused  by  diffraction  from  the  edges  of  the 
output  mirror.  Other  difficulties  are  cuased  by  a  transverse  spatial  de¬ 
pendence  on  gain  and  a  spatial  dependence  on  the  index  of  refraction  of 
the  laser  medium.  Several  previous  attempts  to  develop  an  analytical 
model  have  been  applicable  to  only  specific  problems  or  excessive  in 
computational  costs.  The  geometrical  technique  presented  by  Moore  and 
McCarthy  (Ref  5)  is  very  inexpensive  but  ignores  the  edge  diffraction 
effects  and  is  therefore  not  very  accurate  at  low  Fresnel  numbers.  The 
iterative  technique  presented  by  Rensch  and  Chester  (Ref  7)  does  not 
account  for  transversely  varying  gain,  gives  solutions  for  only  lowest- 
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loss  mode,  and  becomes  extremely  expensive  at  high  Fresnel  numbers.  The 
matrix  technique  suggested  by  Sanderson  and  Streifer  (Ref  8)  also  fails 
to  Include  transverse  gain  and  is  very  expensive  for  higher  Fresnel  numbers 
but  does  allow  for  solutions  corresponding  to  higher  loss  models.  Siegman 
and  Szlklas  (Ref  12)  proposed  a  Fast  Fourier  Transform  technique  which 
Includes  the  effects  of  edge  diffraction,  transversely  varying  gain,  and 
index  of  refraction  buy  only  provides  a  solution  for  the  lowest  loss  mode 
and  despite  its  expense  is  the  widely  used  technique  presently  available. 

Objective 

The  first  objective  is  to  understand  the  technique  developed  by  Moore 
and  McCarthy  (Ref  6)  to  study  positive  branch  confocal  unstable  resonators 
including  a  symmetric  gain  distribution.  This  effort  includes  an  analysis 
of  their  assumptions,  the  limitations  of  their  theory,  and  the  numerical 
integration  techniques  they  employed.  The  second  objective  is  to  then 
adapt  the  computer  program  developed  by  Moore  and  McCarthy  to  the  AFIT 
C DC  6600/CYBER  74  computer  system.  The  results  obtained  will  then  be 
compared  to  Moore  and  McCarthy's  previous  results  to  validate  the  program 
Implementation  at  AFIT.  The  usefulness  of  this  program  as  an  analytical 
tool  will  next  be  demonstrated  by  studying  far  field  distributions  for 
several  different  gain  distributions.  The  final  objective  is  to  extend 
this  program  to  study  the  made  characteristics  for  a  saturated  gain  dis¬ 
tribution.  This  requires  the  calculation  of  the  saturated  gain  distri¬ 
bution  from  the  output  intensity  distribution  by  an  iterative  technique 
to  arrive  at  a  self-consistent  solution. 


Organization 


This  report  contains  five  chapters  and  four  appendices.  The  chapters 
following  this  introduction  are  II,  Theoretical  oackground;  III  Results  of 
Basic  Program;  IV  Analysis  of  Iterated  Gain  Models  and  finally  V  Recomenda- 
tions  and  Conclusions.  A  summary  of  Moore  and  McCarthy's  paper  is  pre¬ 
sented  in  Chapter  II  along  with  the  theoretical  background  necessary  to 
propagate  the  field  to  the  far  field  and  calculate  the  saturated  gain  model. 
The  results  of  using  the  program  in  its  original  form  are  discussed  in 
Chapter  III  including  the  fare  field  calculations.  In  Chapter  IV,  several 
methods  of  deriving  a  saturated  gain  model  are  presented  and  their  vali¬ 
dity  is  discussed.  Chapter  V  contains  specific  conclusions  and  recommenda¬ 
tions  for  future  work  in  this  area.  The  four  appendices  include  derivations 
of  specific  portions  of  Moore  and  McCarthy's  theory  and  explain  the  use  of 
the  computer  program. 
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II.  Theoretical  Background 


This  chapter  contains  a  review  of  the  basic  theory  proposed  by  Moore 
and  McCarthy  (Ref  6)  which  is  used  to  study  the  eigenmodes  of  unstable 
resonators.  This  work  is  hereafter  referred  to  as  Ref  6.  This  discus¬ 
sion  is  designed  to  provide  an  understanding  of  their  theory  when  con¬ 
sidered  in  conjunction  with  Ref  6.  Specific  difficult  sections  of  their 
theory  are  expanded  in  the  appendicies.  Also  included  in  this  chapter 
are  discussions  of  the  method  used  to  find  the  far  field  intensity  dis¬ 
tribution  and  the  relationship  of  the  output  field  to  the  saturated  gain 
distribution. 

Review  of  Moore  and  McCarthy's  Theory 

The  specific  case  that  Moore  and  McCarthy  assumed  was  that  of  a 
positive-branch  confocal  strip  resonator.  The  basic  geometry  of  this 
resonator  is  illustrated  in  Figure  1.  The  resonator  consists  of  a  large 
concave  mirror  and  a  small  convex  feedback  mirror  which  are  aligned  to 
share  a  common  focal  point.  In  this  geometry  the  resonator  magnification 
is  M  =  -  ^  d  ,  where  0  is  the  length  of  the  resonator  and  d  is  the  distance 
from  the  feedback  mirror  to  the  focal  point.  Another  important  parameter 
of  the  resonator  is  the  effective  Fresnel  Number  Fgff  =  ,  where  k 

is  the  wave  number  and  a  is  the  radius  of  the  feedback  mirror.  This 
value  is  often  called  Neq  by  other  authors  and  is  important  because  mode 
separation  shows  a  dependence  on  F^. 

Ignoring  the  diffraction  effects  of  the  feedback  mirrors  results  in 


a  geometrical  solution  for  the  eigenmodes  of  the  cavity.  In  the  limit 


Figure  1.  Basic  Geometry 


of  this  geometrical  intepretation,  the  light  traveling  to  the  left  in  the 
cavity  is  a  cylindrical  wave  emanating  from  the  focal  point  while  the  light 
traveling  to  the  right  is  a  plane  wave.  These  two  types  of  waves  suggest 
the  use  of  a  cartesian  coordinate  system  to  represent  the  right  traveling 
plane  wave  and  a  cylindrical  coordinate  system  centered  at  the  focal  point 
to  represente  the  left-traveling  wave.  When  the  non-geometric  analysis  is 
completed  these  waveforms  retain  their  basic  shape  but  are  modified  by  a 
rapidly  oscillating  amplitude  component  which  is  caused  by  the  diffrac¬ 
tion  effects  of  the  output  mirror.  The  field  inside  the  resonator  is 
therefore  of  the  form 


E  =  f(z,x)  exp 


ik(z+Md) 


g(p.e)  exp 


(1) 


and  obeys  the  scalar  Helmkoltz  equation.  Applying  the  scalar  Helmkoltz 
equation  to  Eq  (1)  and  making  the  paraxial  approximation,  the  following 
equations  for  f(z,x)  and  g(p,e)  are  obtained. 


2ik 

2ik 


9  -  g(p.x)  f(z,x) 

-  g(p.x)  g(p.o) 


+  92  U,sx) 


ax* 


♦  C-2  -  0 

K  99^ 


(2) 

(3) 
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Since  the  total  field  E  must  be  zero  on  the  mirrors,  the  boundary  condi¬ 
tions  for  the  geometrical  solution  to  Eq  (1)  are 


f(o,Mde)  =  (Md)'*5  g(Md,e) 

(4) 

g(d.e)  .(1%f(D.de)e21kD 

,  |  e |  <  a/d 

to 

,  |  e  |  <  a/d 

(5) 

The  geometrical  theory  is  applied  to  a  strip  resonator  by  assuming  that 
the  last  terms  on  the  left-hand  side  of  the  rate  equations,  Eqs  (2)  and 
(3),  are  negligible,  allowing  the  rate  equations  to  be  integrated  directly. 
These  second  derivative  terms  must  be  included  to  account  for  diffraction 
effects  form  the  mirror  edges.  This  necessitates  rewriting  f(z,x)  and 
g(p,e)  in  forms  that  take  into  account  the  cylindrical  diffraction  wave¬ 
lets  emanating  from  the  edges  of  the  feedback  mirror  and  propagating 
through  n  round  trips  of  the  cavity.  The  derivation  of  these  expressions 
for  the  amplitudes  of  f(z,x)  and  g(p,e)  is  contained  in  Appendix  A.  The 
resulting  expressions  are 
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The  amplitudes  fn,  gn,  f,  g,  are  slowly  varying  functions,  white  f  ami 
g  are  rapidly  varying,  expect  ally  at  targe  Since  the  primary  goat 

Is  to  calculate  the  amplitude  of  the  field  at  the  output  mirror,  f(P.\), 
the  values  for  the  amplitude  coefficients  fh(P,\)  and  f(P,s)  must  first 
be  found. 

first  the  expressions  for  f,  fg  (b),  and  g,  tg  (7).  ai'e  substituted 
Into  the  rate  equations,  Iqs  (,')  and  (.1),  to  develop  a  new  set  of  rate 
equations  which  Is  dependent  on  the  direction  of  ray  propagation,  s. 
through  the  cavity 

-■  f„(s,\) 

'^n  "  t»(p.x)  gn(s.\)  ,  etc.  (tO 

where  ''  Indicates  the  directional  derivative  along  the  direction  of 
ray  propagation.  In  addition  to  this  set  of  rate  equations  the  boundary 
conditions  at  the  mirrors  must  be  found  for  each  amplitude  coefficient. 
This  Is  easily  accomplished  for  the  left  end  of  the  cavity  since  It  is 
essentially  an  Infinite  constant  p  surface  over  which  only  one  condition 
must  be  met.  Ibis  condition  Is  that  the  right -travel Ing  wave  must  cancel 
the  left-traveling  wave  on  the  mirror,  lhe  boundary  conditions  at  the 
right  end  of  the  cavity  are  mot'e  difficult  since  two  conditions  must  be 
met.  The  two  conditions  are  that  the  two  fields  must  cancel  each  other 
on  the  mirror  (-a-\-a)  while  g(d,\)  must  be  ?eiv  outside  the  mirror 
(x'|a|).  The  boundary  conditions  are  at  the  left  end  of  the  cavity  aiv 

(on-Md)''1  f„(Nd.xMMd-o  q, (Nd.O  n-l„Vt,..N  (d) 


00) 


f(Md.x)-(Md)‘S  g(Md.x) 

The  boundary  conditions  at  the  right  end  are  then  found  by  using  the 
Horw1t2  approximation  to  the  Fresnel  diffraction  Integral  (Rel  3).  The 
method  used  Is  presented  In  more  detail  In  Appendix  B.  This  numerical 
Integration  yields  not  only  the  boundary  conditions  on  the  right-hand 
mirror  but  an  expression  for  g0(o,x)  which  Is  the  amplitude  of  the 
diffracted  wave  on  Its  first  trip  across  the  cavity.  The  boundary  con¬ 
ditions  at  the  right  end  of  the  cavity  are 

e21kd(t>n-d)'S  fn(d,x)-(d-c  „)'**  g_n(d.x)  n«1 .2.3. . . .N-l  (11) 

et1kd  |  (t'rrd)'  exp  ^  fN(d,x)  +  rN(d. -x )  4  t(d.x)  j 

•  d'*1  g(d,x)  (U1) 

The  equation  for  g0(o.x)  derived  from  the  Integral  Is 

t 

n»T 

x  exp  S  Vk(?3_<l)  (  +  )  - 

■  n  u 

-exp  11  fn(d.-a)(.>n ♦  f(d.a)  '  (13) 

The  right  side  of  this  equation  must  be  multiplied  by  the  exponential  of 
the  gain  Integral  In  order  to  be  valid  over  the  whole  length  of  the  cavity. 
In  this  equation  ;>  Is  the  angle  with  respect  to  the  axis  of  the  ray  from 
the  mirror  edge  (d.a)  to  the  point  where  g0  Is  evaluated. 
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Figure  2.  Illustration  of  3  and  >’n(  +  ) 


The  angles  3n(+)  are  the  angles  with  respect  to  the  axis  of  the  rays 
emanating  from  the  virtual  image  points  (;>_n,  -x_n)  and  striking  the 
output  mirror  edge  at  (d,a)  (i.e.  jj  (lVn)"!). 

The  expressions  for  the  amplitude  coefficients  fn(d,x)  and  f(d,x) 
at  the  output  mirror  are  now  developed  in  order  to  find  the  total  out¬ 
put  field  f(d,x).  The  g0(o,x)  term  is  propagated  back  and  forth  through 
the  cavity  by  alternatively  solving  the  rate  equations,  Eq  (8),  and  apply¬ 
ing  the  boundary  conditions,  Eqs  (9)  and  (11).  After  a  significant  amount 
of  algebra  the  following  equation  for  fn(d.x)  is  obtained 


f„(d,x)  -  e2,("-|>kt’  Sf>- 


f  1  (d  ,Sf|  ( x ) ) 


where  sn(x)  is  a  type  of  shorthand  notation  which  represents  the  path  a 
ray  must  follow  in  order  to  end  up  at  (d.x)  after  making  n  round  trips 
through  the  cavity.  The  expression  for  sn(x)  is  developed  in  Appendix  C. 
The  function  n(x)  is  the  exponential  of  the  integrated  gain  mirror 


(d,a),  makes  n  round  trips  across  the  cavity  then  strikes  the  output 
mirror  plane  at  a  point  (d,x).  It  is  given  by  the  following  expression 


rn 


M]-m 
1-M-  n 


(-M'nx  +  a) 


p'  M 


-1-n+m 


l  K  " _ 

d  1-M-Zn 


{x-M-"a))dp-  *  f">  ( 


+  (-N  nx  +  a))  dp- 


-n+m 

Y  (x-M  a) 
n 


(15) 


An  expression  for  fj(d,x)  is  developed  by  propagating  g0(p,x)  t.o  the 
left  end  of  the  cavity.  Upon  examination  of  the  expression  for  f](d,x) 
the  major  transverse  dependence  is  contained  in  the  term  rl^/a_x/m)  which 
is  factored  out  allowing  all  the  remaining  slowly  varying  transverse  terms 
to  be  collected  in  a  new  term  called  q(x).  The  resulting  expression  for 
f  1  (d ,x )  becomes 


f] (d,x) 


e2lkD  Wg  r]  (x) 


bM. 


(a-x) 

m 


(16) 


The  expression  for  q(x)  containing  all  the  slowly  varying  transverse 
terms  becomes 

(exp  (-^  {£) 

m=l 


X  M* /Z  T  m(a)q  ( 


M'W  'm 

T+FFm 


a) 
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exp 


x 


*( 


/ 1 ka2  1+H“m 


)  «2tmk0  p/Z  r  m(.s) 


+d  a  f(d,a)(l-M"2) 


(in 


In  order  to  solve  this  equation  for  the  modes  of  the  resonator,  an  expres¬ 
sion  for  f(d,a)  must  be  found.  This  expression  Is  found  by  applying  the 
boundary  conditions  at  both  ends  of  the  cavity  and  then  using  the  approxi¬ 
mation  to  the  exponential  terms  for  large  N.  For  large  N  the  exponentials 
1  ka2 

become  exp  (-jj  ).  The  resulting  expression  for  f(d,x)  Is 

f(d,Mx)  •  [  d‘**  exp  mN  |rN(x)  +  r'N(-x)| 

X  q(£)  a'1  (l-M-2)"1  +  f(d,x)  ]  v> 

/  rMd  /'Md  y  •  V 

X  exp(  J  G-  (p*  ,Mx)dp '  +  J  G‘  (p*t  *§-  )  dp*)  (18) 

In  Its  limiting  form  for  large  N,  the  expression  for  r‘N(x)  approaches 
a  symmetric  function  so  the  solution  for  the  antisymmetric  modes  become 

A 

f(d,x)  ■  o  .  For  the  symmetric  modes  the  solution  to  Eq  (16)  becomes 


X  2  T*n(x)  q  (*)  a"1  (1-M‘2f 1  (19) 

Using  this  equation  evaluated  at  xaa  In  Eq  (14)  yields  an  equation  In  which 
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the  only  unknown  amplitude  is  q(x).  The  final  approximation  necessary 
in  order  to  get  a  polynomial  equation  for  y  is  to  assure  that  q(x)  is 
constant  for  all  x.  Moore  and  McCarthy  discuss  two  approximations 
which  could  be  used  to  accomplish  this.  Both  approximations  appear 
to  be  valid  assuming  the  major  transverse  dependence  has  already  been 
separated  out  of  the  expression  for  q(x). 

The  final  polynomial  expression  for  y  is  obtained  by  using  the  above 
symmetric  mode  expression  for  f(d,a)  and  then  dividing  by  q(x)  to  get 
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where  ■»  TT+ftpnf  .  This  is  the  equation  which  is  solved  numerically 
to  calculate  the  roots  of  the  polynomial  corresponding  to  the  eigenvalues 
of  the  various  eigenmodes  of  the  resonator. 

Once  the  roots  of  y  are  found  the  amplitudes  fn(d,x)  and  f(d,x)  can 
be  calculated  and  substituted  into  Eq  (16)  to  calculate  the  field  at  the 
output  mirror.  The  resulting  field  can  only  be  found  out  to  the  shadow 
boundary,  x*Ma,  because  of  the  divergence  which  occurs  at  the  end  points 
of  the  diffraction  Integral  used  to  solve  for  the  boundary  conditions 
at  the  right  hand  mirror. 

Moore  and  McCarthy  chose  to  extend  the  calculation  for  the  field 
through  the  shadow  boundary  by  using  a  slightly  different  form  of  the 
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ntegral  equation  used  to  solve  for  the  right  hand  boundary  conditions. 
This  new  form  of  the  Integral  can  be  expressed  as  a  Rublnowlcz  line 
Integral.  Their  changes  of  variables  and  use  of  the  Rublnowlcz  line 
Integral  have  the  effect  of  pushing  the  discontinuity  outside  of  the 
region  of  Interest.  The  resulting  equation  for  the  output  field  which 
Is  used  by  the  computer  program  is 
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Propagation  to  the  Far  Field 

The  purpose  of  this  section  Is  to  discuss  the  method  which  Is  used 
to  determine  the  far  field  pattern. 

As  outlined  In  the  previous  section,  the  output  field  of  the  reson¬ 
ator  can  be  calculated  and  Is  denoted  f(D,x).  The  output  field  which  Is 
actually  propagated  Is  this  field,  f^(D,x),  minus  the  central  portion 
that  is  reflected  back  into  the  cavity  by  the  output  mirror.  For  the 
purpose  of  this  discussion,  the  field  to  be  propagated  Is  denoted  U(x). 
Since  the  far  field  pattern  Is  desired,  the  Fraunhofer  diffraction  equa¬ 
tion  Is  used  to  propagate  the  field  *(x)  to  the  observation  plane  accord¬ 
ing  to  the  formula  (Ref  1:61) 
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In  this  analysis  the  one-dimensional  Fourier  transform  is  defined  as 

oo 

F|g(x)l  *  J^)e'1MxX  dx  (23) 

where  g(x)  is  a  complex  function  and  fx  Is  an  Independent  variable 
called  a  spatial  frequency.  For  the  case  at  hand,  the  spatial  frequency 
term,  fx,  is  defined  as 
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(24) 


The  last  portion  of  the  Fraunhofer  diffraction  integral  can  r.  ,  be 
written  as 
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dx  « 


The  Fraunhofer  equation  for  the  far  field  then  becomes 
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(26) 


The  complex  far  field  distribution,  U(xQ),  can  therefore  be  found 
directly  by  taking  the  Fourier  transform  of  the  near  field,  U(x). 


Saturated  Gain 

The  magnitude  of  the  laser  threshold  or  steady  state  oscillation 
condition  for  a  stable  resonator  is  defined  as 


1  =  [LOSS]  e2gD  (27) 

where  LOSS  is  the  sum  of  all  the  loss  terms  including  media  and  output 
losses  and  g  Is  the  gain  necessary  to  support  steady  state  oscillations. 
For  an  unstable  resonator  this  equation  becomes 

u  =  |  LOSS  ]  e2gD  (28) 

where  y  is  the  root  of  the  polynomial  expression  which  was  previously 
determined.  This  equation  indicates  the  losses  per  round  trip  must  be 
offset  by  the  gain  in  the  laser  medium  In  order  to  sustain  laser  oscilla¬ 
tions.  The  gain  of  the  medium  is  not  constant  but  depends  on  various 
material  parameters  and  the  field  Intensity  as  follows 
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g 


(29) 


where  gQ  Is  the  small  signal  gain  coefficient  of  the  medium,  Isat  is  the 
intensity  level  necessary  to  cause  the  population  inversion  of  the  laser 
to  drop  to  one-half  its  nonsaturated  value,  and  I  is  the  intensity  in 
watts  per  square  meter  of  the  optical  field  (Ref  13:108).  The  intensity 
can  be  determined  by  squaring  the  magnitude  of  the  output  field,  f(D,x). 
In  the  present  analysis  the  rapid  transverse  variations  of  the  field 
and  therefore  the  intensity  cause  an  expected  transverse  variation  in  the 
gain  so  that  Eq  (29)  becomes 


9(x) 
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WL 
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(30) 


where  (x)  indicates  the  transverse  dimension.  This  equation  is  only  valid 
for  homogeneously  broadened  media,  however,  many  present  applications 
utilize  this  type  of  media.  If  the  desired  medium  is  inhomogeneous ly 
broadened,  a  different  equation  would  have  to  be  utilized. 
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III.  Results  of  Basic  Program 


The  purpose  of  this  chapter  is  to  discuss  the  use  of  the  basic 
computer  program  without  the  changes  necessary  to  include  a  saturated 
gain  profile.  The  changes  necessary  to  make  the  program  compatible  with 
the  AFIT  CDC  6600/CYBER  74  computer  system  are  discussed.  Representative 
plots  of  the  magnitude  and  phase  of  the  output  field  are  presented  and 
compared  to  expected  results.  Finally  the  usefulness  of  this  program  as 
an  analytical  tool  is  demonstrated  by  changing  the  program  to  calculate 
this  for  field  distribution. 

Validation  of  Progt am 

The  program  that  was  provided  by  AFWl  had  undergone  several  signifi¬ 
cant  modifications  since  being  initially  developed  by  Moore  and  McCarthy. 
The  primary  modification  was  that  two  programs  ROOTS  and  MODES  were  com¬ 
bined  into  one  program  thereby  allowing  the  direct  calculation  of  the 
output  field  for  a  selected  number  of  lowest-loss  modes.  The  plotting 
routines  for  both  the  original  program  and  AFWL  modifications  were  not 
compatible  with  the  AFIT  system  and  therefore  required  significant  modi¬ 
fication.  The  ROOTS  program  had  also  been  modified  to  use  a  neh  complex 
root  solving  routine  ZRPCC  instead  of  the  CPQR  routine  which  was  used  by 
Moore  and  McCarthy.  In  addition  the  loaded  resonator  portion  of  the 
program  had  not  been  utilized  since  the  two  programs  were  combined  and 
did  not  calculate  the  expected  field  distributions  when  it  was  first 
utilized. 

The  program  was  methodically  checked  out  to  correct  these  areas. 
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Table  1. 

Validation 

Cases 

Feff 

Mag 

Case 

Previous  Program 

Real  Imag 

Current  Program 

Real  Imag 

8.4 

2.0 

8are  Cavity 
Lowest  loss 

2nd  lowest 

0.95)72 

* 

-.006 02 

* 

0.96175 

1.09313 

-.005947 

-.779319 

Positive  Gain 

0.97160 

-.003881 

0.97076 

-.003932 

16.4 

2.9 

Bare  Cavity 

0.96145 

-.00427 

0.96143 

-.004147 

Positive  Gain 

0.97999 

-.00280 

0.97936 

-.002798 

Negative  Gain 

* 

★ 

0.94903 

-.004514 

*  no  previous  example 


First  all  the  unused  portions,  such  as  old  plotting  routines  and  duplica¬ 
ting  subroutines,  were  deleted.  A  new  plotting  routine,  HGRAPH,  was  added 
which  plots  the  output  field  intensity  and  phase  for  a  specified  number  of 
lowest  loss  modes  of  the  cavity.  The  gain  routines  that  calculated  the 
gaussian  gain  distributions  were  modified  to  produce  similar  results  to 
those  contained  in  Ref  (6).  The  input  and  output  statements  were  also 
modified  to  allow  easy  use  of  the  program  from  the  INTERCOM  terminals. 

The  specific  cases  presented  in  Table  1  were  used  to  facilitate 
comparison  to  the  results  obtained  by  Moore  and  McCarthy.  The  table 
includes  the  complex  values  of  the  root  of  the  polynominal,  u,  from 
both  Moore  and  McCarthy's  paper  (Ref  6}  and  the  program  as  it  is  presently 
configured  on  the  AFIT  computer  system.  There  are  slight  differences  in 
the  values  of  the  calculated  roots.  The  largest  difference  between  cal- 

.4 

culated  real  parts  is  8.4  x  10  which  is  <0.1%  error  while  the  largest 


I 


18 


difference  between  imaginary  parts  is  1.23  x  lO'4  which  is  ;  V*  error. 
These  very  slight  differences  are  attributed  round  off  errors  and  the 
use  of  the  new  complex  root  finding  routine  XRPCC  instead  of  the  slower 
CPQR.  The  larger  error  for  the  imaginary  term  is  less  important  when 
the  magnitude  of  the  root  is  considered. 

The  program  has  also  been  configured  to  provide  plots  of  the  magni¬ 
tude  and  phase  of  the  output  field.  The  plots  for  the  first  case  in 
Table  1  are  included  here  while  the  plots  for  the  rest  of  the  cases  are 
contained  in  Appendix  D.  These  plots  can  be  compared  directly  to  Figures 
13,  11,  8  and  10  from  Ref  6  but  could  not  be  reproduced  well  enough  for 
inclusion.  After  a  qualitative  comparison  of  the  plots  from  Ref  6  and 
the  present  program,  the  curves  are  seen  to  have  the  same  height  and 
variations  in  the  transverse  dimension.  An  additional  check  of  the  gain 
portion  of  the  program  was  obtained  by  plotting  the  calculated  values 
for  gamma  prime  versus  the  transverse  output  dimension.  This  plot 
(Figure  5)  can  be  compared  directly  to  Figure  5  of  Moore  and  McCarthy's 
paper  (Ref  6).  The  plot  shows  that  the  gamma  function  approaches  a 
symmetric  function  for  an  increasing  number  of  passes  which  corresponds 
to  the  case  of  large  N.  These  Independent  checks  indicate  that  the 
program  is  correctly  calculating  and  plotting  the  output  fields  for 
both  bare  and  loaded  cavities. 

Far  Field  Patterns 

The  program  was  next  modified  to  calculate  the  magnitude  of  the 
far  field  distribution  by  using  the  Fourier  transform  technique  discus¬ 
sed  in  Chapter  II.  The  Fourier  transform  is  calculated  by  a  complex 
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Fast  Fourier  transform  routine,  FFT2C,  from  the  IMSL  program  library. 
Since  the  basic  program  only  calculates  half  of  the  output  field,  the 
field  must  be  mirrored  around  the  optical  axis  to  represent  the  other 
half  of  the  field  and  then  have  the  center  portion  deleted  to  represent 
the  obscuration  introduced  by  the  output  mirror.  The  shape  of  the  output 
field  which  is  finally  propagated  is  illustrated  in  Figure  6.  The  trans¬ 
verse  dimension  in  the  near  field  is  x/a  so  that  the  mirror  edge  is 
normalized  to  1.  The  transverse  dimension  in  the  far  field  becomes 


by  the  scaling  property  of  the  Fourier  transform  when  the  near  field  is 
normalized  to  the  radius,  a,  of  the  small  mirror. 

The  near  and  far  field  intensity  distributions  for  the  two  lowest 
loss  modes  of  a  loaded  resonator  with  F  ^  *  8.4  and  a  magnification 
*  2.0  are  presented  in  Figures  7,  8,  9,  and  10.  It  is  obvious  that 
both  fare  field  patterns  have  lost  any  remanent  of  the  centered  obscur¬ 
ation  and  are  focused  down  to  a  relatively  small  area.  The  sharp  dis¬ 
continuities  caused  by  the  centered  obscuration  tend  to  make  the  output 
field  approach  a  square  wave  with  relatively  small  variations  in  ampli¬ 
tude  until  the  shadow  boundary  is  reached.  When  the  limiting  case  of 
a  pure  square  wave  Is  propagated  the  far  field  distribution  becomes  a 
delta  function.  This  property  of  a  square  wave  coupled  with  the  general 
shape  of  the  output  field  accounts  for  the  relatively  minor  differences 
between  the  lowest  loss  and  next  lowest  loss  far  field  patterns.  A 
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somewhat  more  detailed  study  of  the  far  fields  could  be  accomplished  by 
calculating  the  Intensity  within  the  first  Airy  Disk.  This  Is  often 
called  the  "power  In  a  bucket"  and  allows  a  more  quantitative  study  of 
the  far  field  distribution  than  simply  observing  the  far  field.  Regard¬ 
less  of  the  method  used  to  study  the  far  field  patterns,  the  utility  of 
the  basic  program  was  demonstrated  by  allowing  quick  calculations  for 
several  different  resonator  parameters. 
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IV.  Analysis  of  Iterated  Gain 

The  purpose  of  this  chapter  is  to  discuss  the  modifications  to  the 
basic  program  necessary  to  calculate  a  saturated  gain  profile.  The 
application  of  the  theory  to  the  parameters  calculated  by  the  program 
will  be  discussed  first  followed  by  a  review  of  several  attempts  at 
developing  reasonable  models.  This  Chapter  is  concluded  with  a  discus¬ 
sion  of  possible  problem  areas  which  must  be  addressed  in  order  to 
derive  an  ideal  saturated  gain  model. 

Applications  to  Present  Program 

Moore  and  McCarthy's  basic  technique  can  calculate  the  eigenmodes 
for  either  a  "bare"  cavity  or  a  cavity  which  has  an  arbitrary  symmetric 
gain  distribution.  They  chose  to  provide  examples  where  the  slowly 
varying  gain  function  was  either  a  positive  or  negative  Gaussia  dis¬ 
tribution.  However,  in  a  real  laser  resonator  the  gain  distribution 
is  neither  constant  nor  arbitrary  but  is  actually  dependent  on  the 
field  intensity  according  to  the  relationship  developed  in  Chapter  II. 
This  relationship  of  gain  to  Intensity  is  repeated  here 


G(x) 
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(32) 


where  gQ  is  the  small  signal  gain  of  the  laser  medium,  (lx)  is  the  out¬ 
put  Intensity  at  a  distance  x  from  the  optical  axis,  and  Isat  is  the 
saturation  Intensity  which  Is  dependent  on  pumping  levels  and  mode 
volume. 
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Figure  11.  Iteration  Flow  Diagram 


To  calculate  the  required  gain  distribution  the  following  iterative 
approach  illustrated  in  Figure  11  was  used  to  arrive  at  a  self  consistent 
solution  for  the  output  field  and  the  gain  distribution.  A  bare  cavity 
was  assumed  on  the  first  iteration  to  determine  an  initial  intensity 
distribution.  A  new  gain  distribution  was  then  calculated  and  used  in 
turn  to  find  new  values  for  y  the  root  of  polynomial.  This  new  root 
is  then  used  to  calculate  a  new  output  field  from  which  a  second  gain 
distribution  was  calculated.  This  process  was  repeated  until  the  value 
for  y,  the  gain  distribution,  and  the  output  field  distribution  changed 
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by  less  than  0.1%  from  one  iteration  to  the  next.  Convergence  using  this 
criteria  always  occurred  within  10  iterations  and  usually  within  5  itera 
tions. 

The  gain  distribution,  G(x),  calculated  using  Eq  (22)  could  not  be 
used  directly  since  the  program  had  been  arranged  to  use  the  related 
value  G’(x)  which  is  defined 

G*(x)  =  G(x)  -  Q/D  (33) 


where 


Q 


,o)  dp* 


(34) 


or  specifically  for  this  case 

Q  =  G(o)D  (35) 

The  formula  used  to  calculate  a  new  G'(x)  is  therefore  taken  to  be 

G*(x)  =  G(x)  -  G(o)  (36) 


or 
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By  Inspection  it  is  easy  to  see  that  G'(o)  =0  as  is  required  by  its 
definition. 

This  approach  of  using  only  the  field  at  the  output  mirror  does  not 
take  into  account  the  dynamic  changes  that  occur  during  a  round  trip  of 
the  cavity.  The  approach  is  useful  however  in  determining  the  validity 
of  calculating  a  saturated  gain  distribution  using  this  iterative  ap¬ 
proach  is  useful  however  in  determining  the  validity  of  calculating  a 
saturated  gain  distribution  using  this  iterative  approach.  This  method 
could  be  extended  to  calculate  the  intensity  and  gain  distribution  at 
several  locations  throughout  the  cavity  to  account  for  the  dynamic 
changes  in  the  gain. 

Attempted  Models 

In  the  above  expression  for  G'(x),  the  value  for  Isat  is  a  material 
dependent  parameter  that  cannot  be  fixed  without  assuming  a  specific 
lasing  material  and  pumping  level.  Since  the  goal  is  to  achieve  a  model 
for  saturated  gain  without  having  to  specify  the  lasing  material  and 
since  Isat  is  simply  a  constant  that  affects  the  gain  uniformly  it  was 
at  first  ignored.  In  this  case  the  equation  for  G'(x)  became 

G  ^  =  9o(  1  +  I(x)  "  1  +  !(o)  )  (38) 

The  value  of  gQ  is  also  dependent  on  the  media,  but  since  the  program 
uses  round-trip  gain  coefficients,  it  is  considered  a  dimensionless 
quantity  representing  the  gain  per  pass.  It  was  chosen  to  be  a  constant 
between  0  and  10  where  0  corresponds  to  below  threshold  conditions  and 
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10  corresponds  to  very  large  gain  which  should  cause  the  output  field 
to  saturate. 

After  several  runs  using  this  technique.  It  became  obvious  that  the 
large  values  of  output  intensity  (~  1000)  completely  dominated  the  expres¬ 
sion  for  G ‘ (x )  and  closely  approached  the  bare  cavity  (G * { x )  -  0)  solu¬ 
tion.  This  large  Intensity  also  completely  negated  any  variation  In  the 
gQ  term.  It  was  therefore  obvious  that  some  value  for  Isat  had  to  be 
chosen  In  order  to  calculate  a  reasonable  value  for  saturated  gain.  To 
avoid  choosing  a  media,  Isat  was  chosen  to  be  simply  the  Intensity  along 
the  optical  axis  of  the  resonator.  The  equation  for  G  (x)  now  became 


G'(x) 


(39) 


with  gQ  rental ng  as  before.  This  method  allowed  the  calculation  of  a 
gain  distribution  with  a  reasonable  variation  but  since  I(o)  was  dif¬ 
ferent  for  each  Iteration  the  assumed  value  for  Isat  was  not  very 
realistic.  In  addition,  this  method  did  not  exhibit  the  desired  decay¬ 
ing  modes  for  gain  less  than  threshold  or  saturated  Intensity  for  ex¬ 
tremely  large  gain  coefficients. 

This  method  was  modified  to  use  the  first  Iteration  value  of  I (o ) 
as  the  value  of  I$at  all  future  Iterations.  This  Is  also  an  artifi¬ 
cial  way  of  selecting  Isat  and  still  has  the  problem  of  not  acting 
properly  for  varying  gain. 

For  the  final  attempt,  the  gain  was  fixed  at  the  value  required 
for  steady-state  oscillations  according  to  the  formula 
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u  »  M-*5  e2gD 

Solving  this  equation,  the  gain  necessary  to  support  steady-state 
oscillations  Is  seen  to  be 


(40) 


g$  =  .laJn-M..+.  \  l.n  M 


(41) 


where  u  Is  the  magnitude  of  the  root  of  the  lowest  loss  mode  and  D  is 
the  length  of  the  cavity.  The  value  for  Isat  was  then  taken  to  be 
which  assumes  that  Isat  Is  some  fraction  of  the  Intensity  along  the 
optical  axis.  This  Is  again  an  artificial  value  because  the  laser  media 
Is  not  considered,  but  does  account  for  increased  pumping  by  increasing 
the  value  for  n.  The  equation  for  G  (x)  in  this  case  becomes 


For  small  values  of  n  (less  than  5)  this  method  caused  convergence  but 
still  did  not  correctly  account  for  other  than  threshold  gain  assumptions. 
For  larter  values  of  n,  the  n  term  dominated  the  equation  and  the  gain 
distribution  apporached  the  bare  cavity  level. 

Problem  Analysis 

None  of  these  attempts  at  calculating  realistic  naturated  gain  dis¬ 
tributions  were  adequate  in  all  respects,  but  some  Important  conclusions 
can  be  derived.  First,  a  specific  laser  medium  must  be  specified  in 
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order  to  derive  the  expressions  for  the  small  signal  gain,  gQ,  and  the 
saturation  Intensity,  ISat-  Without  specifying  the  medium  these  values 
are  set  rather  arbitrarily  to  be  of  the  proper  magnitude.  The  second 
and  primary  conclusion  is  that  the  Iterative  routine  does  calculate  an 
appropriate  shape  for  the  transverse  gain  distribution  but  does  not 
account  for  saturation  effects  when  1 t  recalculates  the  output  field  dis- 
trlbutlon.  The  basic  problem  Is  that  the  ROOTS  subroutine  always  must 
assume  a  steady-state  condition  in  order  to  calculate  the  value  of  u. 

When  the  program  goes  on  to  calculate  the  output  field  a  steady-state 
condition  has  already  been  assumed  so  the  expected  decay  for  gain  lower 
than  the  threshold  value  does  not  cause  the  output  field  to  decay  as 
would  normally  be  expected.  The  proposed  solution  is  to  use  a  different 
value  for  u  when  the  output  field  is  calculated  that  the  value  used  to 
find  the  roots  of  the  polynomial.  One  possible  method  for  accomplishing 
this  "decoupling"  is  outlined  below  but  has  not  been  implemented  in  the 
program.  As  the  program  is  presently  configured,  the  magnitude  of  the 
steady  state  gain  is  defined  by  the  equation 

u  -  M'VG(o)L  (43) 

on  the  first  iteration  through  the  program.  To  account  for  the  satura¬ 
tion  effects  on  subsequent  iterations  through  the  program,  a  new  value 
of  vi  must  be  defined  by  the  equation 
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Solving  for  u$  in  terms  of  y,  which  Is  the  basic  root  of  the  polynomial, 
the  following  expression  is  obtained 


This  value  can  now  be  used  In  Eq  (20)  to  calculate  the  output  field  which 
should  now  be  decoupled  from  the  steady  state  assumption  necessary  to 
calculate  the  first  u.  This  procedure  should  be  validated  by  changing 
the  program  then  applying  either  very  low  or  high  gain  factors  to  observe 
the  effect  on  the  shape  of  the  output  field. 
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V.  Conclusions  and  Recommendations 


Conclusions 

The  first  objective  was  accomplished.  The  theory  presented  by 
Moore  and  McCarthy  was  reproduced  in  detail  and  the  necessary  approxi¬ 
mations  were  analyzed.  The  Appendicies  A,  B,  and  C  and  be  used  to  under¬ 
stand  some  of  the  particularly  difficult  sections  of  Moore  and  McCarthy's 
theory  without  having  to  go  back  to  basic  principles.  The  second  objec¬ 
tive  was  also  met.  The  program  is  now  completely  operational  locally  and 
has  been  modified  to  do  some  basic  resonator  analysis.  The  program  can 
now  be  used  to  study  both  loaded  and  bare  cavities  for  a  variety  of 
effective  Fresnel  numbers  and  magnifications  with  very  low  cost. 

The  third  objective  was  not  fully  achieved,  however,  significant 
progress  was  made  toward  isolating  the  problem  area.  As  presently  con¬ 
figured  the  basic  shape  of  the  saturated  gain  distribution  can  be  found 
and  a  method  for  accounting  for  saturation  effects  has  been  identified 
but  not  validated.  Once  this  problem  is  overcome  the  mode  characteristics 
of  resonators  with  realistic  saturated  gain  distribution  can  be  accom¬ 
plished  at  realtively  low  cost. 

Recommendations 

The  analysis  of  the  application  of  saturation  effects  should  be 
continued  to  derive  an  acceptable  technique  that  behaves  properly  for 
all  gain  levels.  Once  this  method  is  complete  a  model  for  the  longi¬ 
tudinal  variation  in  gain  should  be  included  rather  than  the  simple 
single  gain  sheet  model  discussed  in  this  paper. 
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The  study  of  the  far  field  distributions  should  be  expanded  to 
allow  "power  In  the  bucket"  calculations.  This  would  be  particularly 
Interesting  for  field  distributions  calculated  using  an  ideal  saturated 
gain  model.  Another  interesting  exercise  would  be  to  expand  this  theory 
and  associated  program  to  a  three-dimentlonal  cavity.  This  would  more 
closely  represent  a  "real"  laser  cavity  and  would  permit  variations  of 
gain  In  two  dimensions  to  simulate  a  flowing  medium. 
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Appendix  A 

Development  of  Amplitude  Equations 


This  appendix  will  cover  the  approximations  necessary  to  arrive  at 
the  expressions  for  the  amplitudes  f  and  g  in  Eqs  (20)  and  (21)  of  Moore 
and  McCarthy's  paper  (Ref  6). 

The  geometrical  theory  predicts  a  plane  wave  for  the  right-traveling 
portion  fo  the  field  inside  the  cavity  and  a  cylindrical  wave  emanating 
from  the  focal  point  for  the  left-traveling  wave.  When  edge  diffraction 
effects  are  included  in  the  right-traveling  wave  is  found  to  be  the  sum¬ 
mation  of  the  basis  plane  wave  and  cylindrical  wavelets  emanating  from 
the  vitual  Images  at  the  points  specifiec  by  (pn,i  xn)  where 


Pn  *  dM2n 

(A-l) 

xn  =  aMn 

(A-2) 

as  specified  in  Ref  6.  The  objective  of  the  following  discussion  is  to 
derive  the  approximations  necessary  in  the  magnitude  and  phase  of  these 
cylindrical  wavelets. 

The  expression  for  a  cylindrical  wave  is  given  by  the  following 
formula  (Ref  2:25) 


♦(r.x)  «  -J-  A(P,x)e1kR 

/IT  (A-3) 

where  A(p,x)  Is  the  amplitude  of  the  wavelet  and  R  is  the  distance  from 
the  virtual  image  to  the  observation  point  as  shown  in  Figure  A-l. 
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Figure  A-l.  Geometry  for  right  traveling  wave 


In  the  magnitude  portion  of  Eq  (A-3)  the  distance  R  may  be  approximated 

R  =  zn  -  X 

:  Pn  -  p  (A-4) 

since  the  paraxial  approximation  has  already  been  assumed.  The  approxi¬ 
mation  for  the  distance  R  in  the  phase  term  requires  more  care.  From 
geometry  it  is  obvious  that 

R  •  -  0)2  *  (x-*/  ]  S  (A-5) 

Using  the  binomial  series  expansion  for  the  square  root 


R 


+  .  .  .J  (A-6) 


The  Fresnel  degree  of  approximation  assumes  that  only  the  first  two  terms 
of  the  series  must  be  used,  therefore. 


F  ■  «<*,  -  ■>  ♦  r  ^ 


(A-7) 


Now  using  the  above  approximations  for  R  in  Eq  A-3  the  field  at  the 
observation  point  due  to  the  n*^  virtual  image  becomes 

.  v  /ik  a-x,)1) 

Lk(fr-d  7T? ) 


(p ,x)  =  (P  -p)-*5  A  (p,x)  e 


(A-8) 


Since  the  exp(ik(pn  -p)  term  is  slowly  varing  it  is  lumped  together  with 
the  amplitude  term  to  give  a  new  slowly  varying  amplitude  expression 

fn(p.x)  «  An(p,x)  eiMr'-f) 

The  field  due  to  the  nth  source  is  therefore 

/i*  l*-*ii)2 

♦n(p»x)  *  (Pn  -p)"1*  fn(p.x)  'eZ  pn  _p 

To  arrive  at  the  total  expression  for  the  amplitude  of  the  right-travel¬ 
ing  wave  the  effect  of  the  fields  from  N  virtual  images  must  be  summed 
along  with  the  basic  plane  wave  to  give  the  following  expression 


(A-9) 


(A-10) 
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-exp  (H  ^(>!^")?)  yp.-«)]  ♦  f(p.x)  (A-ll) 

A 

where  f(p,x)  is  the  amplitude  due  to  the  basic  plane  wave.  This  corres¬ 
ponds  to  Eq  (20)  of  Ref  6. 

This  procedure  is  repeated  to  find  the  amplitudes  for  the  left 
traveling  wave.  The  virtual  images  are  now  to  the  right  of  the  cavity 
as  illustrated  in  Figure  A-2.  The  approximation  for  the  distance  R  is 
found  in  the  following  manner.  From  geometry 


R2  -  (*-*/  *  <x-x/ 


(A-12) 


r 


r 


using  the  fact  that  z-  p2-x2  and  z  =  y  2  -  x  2  and  squaring  all  terms 

n  n  n 


R2  3  P2  +  p2  2  2(p2-x2)  -2xx„ 


substitute  x=ep  and  *n=0nPn  in  the  previous  equation  to  get 


R2  -  p2  ♦  »2  -  2eenppn  -  2ppn  ■'()-  2)(t-  2) 


The  value  under  the  square  root  becomes 


■OTT?)  -  77n 


(A-13) 


(A-14) 


(A-15) 
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Figure  A-2.  Geometry  for  left  traveling  wave 

Since  the  terms  on  the  order  of  0  are  much  smaller  than  the  remaining 
terms,  this  equation  is  approximated  as 


:  i 


+  e2„ 


(A-16) 


when  the  Fresnel  degree  of  approximation  is  assumed.  The  equation  for 
2 

R  now  becomes 


R2  =  (p~Pn)2  +  Ppn(9"en)2 


( A- 17) 


Again  use  the  binomial  series  expansion  for  the  square  root  to  the  Fresnel 
degree  of  approximation  to  get  the  following  expression  for  the  distance  R 


(A-18) 


This  expression  Is  then  substituted  into  the  phase  term  of  a  cylindrical 
wave  Eq  (A-3)  to  yield 
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(A-19) 


r 


>■^1 1 — ■ 


/Ik  (O-Qn)f  ) 

g(p.e)  «  (p-p^’*5  An(P,x)eik(p"pn)  V  rK^T/p/ 

Again  the  exp(ik(p-pn))  term  is  lumped  together  with  the  slowly  varying 
amplitude  coefficient  to  give  a  new  slowly  varying  amplitude  expression. 

gn(p,x)  -  An(P,x)  eik(p’pn)  (A-20) 


The  entire  expression  for  the  left-traveling  amplitude  expression  must 
sum  the  contributions  of  the  main  cylindrical  wave  and  the  wavelets 
emanating  from  all  the  vitual  images  to  the  right  of  the  output  mirror. 
The  final  expression  becomes 


-N+l 

g(p,e) 

n-o )  I )  •  • 


x  gn(p,x)  +  exp 


+ 


g(p»x) 


(A-21) 
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Appendix  B 

Derivation  of  Boundary  Conditions 

The  purpose  of  this  section  is  to  explain  the  procedure  for  solving 
the  integral  equation  and  then  use  its  solution  to  derive  the  boundary 
conditions  at  the  right  end  of  the  cavity  as  well  as  the  expression  for 
the  diffracted  wave  on  its  first  trip  to  the  right. 

As  was  explained  in  Ref  6,  the  equation  which  must  be  solved  to 
match  the  boundary  conditions  at  the  right  end  of  the  cavity  is 


9(p,e)  - 


1  7  2ikn  Is  ad  2ikD 


/  2ikn  V 

2ir  \  1/p  -  1/d  / 


,/f(0,an) 


dn 


Remember  that  from  Eq  (6) 
N 


f(D,an)  * 


■S  M^) 


X  fn(d,-an)  +  f(d,an) 


(8-1 ) 


(B-2) 


Substituting  this  into  the 

9(p-e)  =  k  (wT7d 


integral  expression,  Eq  (B-l)  becomes 
)  ^  ad  e2ikD  (I,+I2+I3) 


(B-3) 
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where 


V  /'I  i  / ilk(an  +x  )2\ 


X  fn(d.a.,)  exp  I  17p-T7d 


(>-  , 
V  i/e-  nr)  *" 


V  ,1  ,  /^Man+xn)2\ 

(en-d)  1  exp  ^  en-T  j 

k  /„  an\2’ 


(  h{0~  ~6Y  \ 

X  fn(d.an)  exp  y  T/;r_-  yd~J  dn 


*1 


z'1  -  /  a'^2\ 

/  f(d.an)  exp  ^  -yr-.-y/y] 


dn 


(B-4) 


(B-5) 


(B-6) 


These  Integral  expressions  can  be  solved  using  an  asymptotic  approxima¬ 
tion  developed  by  Horwltz  (Ref  3).  This  formula  can  be  written  as 
+1 


-T 

Mrf  F(S)  «p  1 

• 

? 

l<(‘-  D\ 

+F(-1)  (  21 (A+B)  )  _1 

exp  (l (A+28+C)) 

+F(1)  (2KA-B))*1 

exp  ( i ( A-2B+C )) 

(B-7) 


For  this  approximation  to  be  valid.  It  Is  assumed  that  A  and  B  are  large 
and  that  B/A  lies  within  the  domain  of  integration  and  not  close  to  n*-l. 


47 


The  function  F(n)  must  also  be  slowly  varying. 
Now  solving  for  using  this  approximation 


•1 


N, 


X  exp  (  i (A^n  -  ?B^n  ♦  Cj))  dn 


(0-0) 


where 


A1 

kac  , 
2d  \ 

(0-9) 

"l 

ka  / 

f  l 

(11-10) 

C1 

■f  ( 

Xn.2k-<L./.  rfli'n.-!*)  ) 

(on-o)(r-o)  ' 

(0-11) 

Solving 

for  the 

partial  tonus  of  the  approximation 

d 

m 

a 

(8-1.') 

C1 

-  »? 
A, 

k 

’  2 

1  ilT,.-3T  (>‘l"  “  t^-vO  | 

(8-18) 

A1 

•  B. 

•  Y 

1  «f*l\v  f  \\L*  ♦  «»d.)  1 

1  rn-d  d(^d)  1 

(8-14) 

A1 

-  ni 

■  Y 

1  P„-d  4o-c0  1 

(8-1 S) 
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A1  ♦  2Bj  +  Cj  ■  | 


1^8} 


!  *  srfl 


( B- 1 6 ) 


A,  -  2B,  ♦  C 


'1 


1  T  L1  7 


l&Sf  *  sferr 


p(a-ed)* 


(B-l  7 ) 


Substituting  these  into  the  approximation  the  solution  for  1^  becomes 
N 


>,  -  I (•„-<>*(&  fc^-ar) 

n«l  n 


(B-l 8) 


Similar  expressions  for  and  1^  can  be  found  yb  using  a  similar  procedure. 
After  finding  expressions  for  and  1^,  Eq  (B-l)  becomes 
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Now  remember  that  the  specific  form  of  the  gi  ,  )  term  on  the  left  side 
of  the  equation  Is  given  by  Eq  (7)  with  n=-m 


x  2, 


g(fd  p’s  (p-p.n)_l5  exp/^-  (8  - 


n=l  ,2 


V 


l/p_n"l/p 


X  g  n(p ,x)  -  exp  / jp  /e  +  x-n 

(  W>-'- 


x)  +  g(p,x)  (B-20) 


After  some  algebra  and  in  the  limit  p-Kl  the  following  exponential  from 
the  left  side 


exp^  ^ 


8  -  _;n 
p-n> 

r/-P~-T/p 


can  be  shown  to  equal  the  following  exponential  from  the  integral  expres¬ 
sion. 


These  exponentials  both  reduce  to 


Hr 


In  the  limited  p-*d.  Similarly,  the  remaining  exponential  on  the  left 
side  can  be  equated  to  an  Integral  exponential  In  the  limit  p-*d  as 
follows 


Now  equating  terms  from  both  sides  that  share  an  equivalent  exponential 
expression,  the  following  partial  boundary  results 


*..-2 


h  /  _  /  ..x  _  1  2irika  *  . 

p  (p-p_n)  9-n(p’x)  ‘  57  dTf/pT/d)"  e 


2ikD 


'pn  l  Va?  d(p-d)  +  p(p_-<d) 


xfn(d-d  dn(i^te;M))  •n'1’2’ . (B‘22) 


In  the  limit  p-kI,  this  equation  becomes 

e21kd(pn-d)"S  fn(d,x)=(d-p_n)"^  gn(d,x)  n=l ,2 . N-l  (B-23) 

This  equation  can  be  duplicated  by  using  the  other  set  of  similar  exponen¬ 
tials.  This  equation  is  the  same  as  Eq  (29)  of  Ref  6.  This  equation  only 
accounts  for  the  terms  from  n*l  through  N-l.  The  Nth  terms  from  the  right 

A 

side  of  the  equation  are  equated  to  the  remaining  g(p,x)  term  on  the  left 
to  give 
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9  (9.x)  ■  ( 


1  f2irik 


■dipN-dl,(p.-dX_ 
ka7  d(p-d)  +  pfp^-c 


.d  e21kD 


fN(d,x) 


X  exp 


ik 

p(xw  -  ed)2 

T 

d(p-a)  +  pIpn- 

ik 

p(xn  +  ed)2 

T 

-  fN(d,-x) 


+  |  ffl  ]  f(d,x)  (B-24) 


In  the  limit  p-+d  and  N  large  this  becomes 

g(d,x)  =  dV1kD  (pN“d ) -Jfi  expj1^2) 

X  (fN(d,x)  t  fN(d,-x)|  +  f(d,x)  (B-25) 

This  is  equivalent  to  the  boundary  condition  Eq  (30)  of  Ref  6. 

The  remaining  terms  are  now  combined  to  give  an  expression  for 
90(p»x)  which  Is  interpreted  as  the  amplitude  of  the  diffracted  wave  on 
its  first  trip  across  the  cavity.  The  expression  is 


These  terms  can  now  be  equated  by  the  following  equivalent  exponentials 


"tlrl-  -IP  »l 


which  are  equated  by  remembering  that  xQ  *  a  and  pQ  *  d. 


(B-28) 


p's(p-d)'‘5  gQ(p,x)  *  ^  ( 
N 

X  £  <vd>'s 


1  /  2irika 


*(vd--*>  I ,ka (s^9 ♦  Sjfiaf1-)  f  ««p|r 

*  fn<d*a)  [«■  Hr 


1  [  2wi ka 


2  ir 


Is  e2,kD  «*+)  |1I-  f 


(B-29) 


This  equation  is  aimilar  to  Eq  (31)  of  Ref  6  and  can  be  made  to  equal  it 
exactly  if  the  variables  3  and  3  ^  are  derived  from  the  coefficients 

A 

of  the  amplitude  terms.  From  the  f(d,a)  term,  the  following  relationship 
must  be  solved  to  find  a  value  for  3. 


(B-30) 


Solving  for  3 


(8-31 ) 


in  -JMVffUll 


Now  using  this  value  in  the  expression  for  the  coefficient  of  the 
fn(d,a)  term 


O, 


<♦> 


-  3) 


a-Xr 


a-ed) 

!p-d) 


a 

1 


-  3 


(B-32) 


Therefore, 


3 


(+)  . 

n 


(B-33) 


Using  the  substitutions  for  xR  and 


»,w 


a 

7 


(1  +  M"n) 


-1 


(B-34) 
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Appendix  C 
Derivation  of  Sn(x) 


The  purpose  of  this  section  is  to  illustrate  the  origin  of  the 
Sn (x )  term.  The  relationship  for  ( x )  first  occurs  in  the  equation 

relating  the  amplitude  fn(d,x)  on  the  nth  pass  to  the  amplitude 
f^(d,S  (x))  on  the  first  pass.  This  equation  is 

(n-D 

f(d.x)  *  JHn-VU  „  ^ 


X 


rfn(x) 
r-,  (S„(x) ) 


f1(d,Sn(x)) 


(C-l) 


In  this  equation  the  Sn (x )  is  a  shorthand  expression  of  the  path  a  ray 
must  travel  in  order  to  end  up  at  x  on  the  final  pass. 

The  procedure  for  deriving  the  general  expression  for  Sp(x)  will 
be  to  first  derive  the  relationship  between  fn(d,x)  and  f  (d,rn(x) ) 
and  then  use  this  to  propagate  all  the  way  back  to  f^(d,sn(x)).  The 
amplitude  f n ( d . x )  will  first  be  propagated  back  to  the  left  hand  mirror 
where  the  appropriate  boundary  conditions  will  be  applied.  This  expres¬ 
sion  will  then  be  propagated  back  to  the  output  mirror  where  the  boundary 
condition  will  again  be  applied.  For  the  purpose  of  this  discussion 
the  exponential  gain  is  assumed  to  be  unity. 

Starting  with  Eq  (32)  of  Ref  6  and  propagating  it  backwards  from 
right  to  left  and  letting  z=D 

fn(d,x)  »  fn(Hd,Xj)  n-1,2 . N  (C-2) 
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where 


Xj  =  x 


M-M 


2n 


1-M 


2n 


-  a 


Mn(M-l) 


n=l  ,2,. . .  ,N 


(C-3) 


Using  the  boundary  condition  of  Eq  (24)  of  Ref  6  at  the  left  end  of  the 
cavity 


fn(Md,V  =(^TTiiT)S  9,.n  (Md. 


Xr) 


I'  n=l ,2,. . . ,N  (C-4) 


The  9-|_n(Md,Xj)  term  is  then  propagated  back  to  the  right  end  of  the 
cavity  using  Eq  (33)  of  Ref  6 

(Sr^T Y  *  ((Sw)  v 


X  gln(d,rn(x))  n=2,3,...,N 


(C-5) 


where 


r„(x) 


=  X 


1-M 


I2  ( 1 '" )  .  . 

a 


1  H-H^ 


n=2,3, . . . ,N 


(C-6) 


Next,  the  boundary  condition  at  the  right  end  of  the  cavity  is  applied 
resulting  in 
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(M2n  M  \  h 

9, .„(«•'„<*»  -  e2,k0 

A  M2(l-n)  M2n_M  \  li 
X  yi^n"^-1  fn-l^d,rn^x^ 

n=2,3 . N  (C-7) 

Using  Eqs  (C-3)  and  (C-6)  the  rn(x)  value  can  be  expressed  in  terms  of  x 
r  (x)  =  x  -  a 

fTn-l  M^n  -1  n=2,3,...,N  (C-8) 

The  relationship  for  fp(d,x)  to  fp  -j (d,rn(x) )  finally  becomes 

fn(d,x)  =  e2ikD  MS  fnl(d,rn(x))  n=2,3 . N  (C-9) 


By  using  a  similar  procedure  the  fn_](d»r  (x))  term  can  be  related  to 
the  fn_2^d,rn-l^x^^  term  t>ie  f°llowin9  manner 


f„.i<d-rn<x»  ■  w-vi'*’” 


This  is  repeated  until 


f2(d,r3(r4(...rn_i(rn(x)  )  ...  )  )  =  e 

x  MS  f^d.S^x)) 


2ikD 


(C-10) 


(C-ll) 
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where 


I 


sn(x)  =  r2^r3^*,rn-l^rn^x^  )•••))  (C-12) 

This  resulting  expression  for  Sn ( x )  is  clearly  the  required  starting 
location  for  a  ray  that  ends  up  at  x  after  n  round  trips. 
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Appendix  D 
Plots 

This  section  contains  the  graphs  for  the  cases  In  Table  1. 
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Figure  D-2.  Fhase,  positive  gain,  lowest  loss  mode 
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Figure  D-3.  Intensity,  bare  cavity , second  lowest  loss  mode 
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Figure  D-5.  Intensity,  bare  cavity , lowest  loss  mode 
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Figure  0-8.  Phase,  positive  gain,  lowest  loss  mode 
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